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Abstract 

We show that the Dirac equation in 3+1 dimensions gives rise to supersymmetric patterns when 
the scalar and vector potentials are (i) Coulombic with arbitrary strengths or (ii) when their sum 
or difference is a constant, leading to relativistic pseudospin and spin symmetries. The conserved 
quantities and the common intertwining relation responsible for such patterns are discussed. 

PACS numbers: 24.10.Jv, 11.30.Pb, 21.60.Cs, 24.80. +y 



1 



The Dirac equation for spin 1/2 particles plays a central role in the relativistic descrip- 
tion of atoms, nuclei and hadrons. In atoms the relevant potentials felt by the electron (or 
muon in muonic atoms) are Coulombic vector potentials. Relativistic mean fields in nuclei 



generated by meson exchanges jl | , and quark confinement in hadrons 2| necessitate a mix- 
ture of Lorentz vector and scalar potentials. Recently symmetries of Dirac Hamiltonians 
with such Lorentz structure have been shown to be relevant for explaining the observed 



degeneracies of certain shell- model orbitals in nuclei ("pseudospin doublets") 



3|, and the 



absence of quark spin-orbit splitting ("spin doublets") |4(, as observed in heavy- light quark 
mesons. The goal of the current letter is to show that the degeneracy patterns and relations 
between wave functions implied by such relativistic symmetries resemble the patterns found 
in super symmetric schemes. The underlying mechanism responsible for such properties will 
be examined. The feasibility of such a proposal gains support from the fact that Dirac 
Hamiltonians with selected external fields are known to be supersymmetric 5|, e.g., for a 
vector Coulomb potential 6|. 

Supersymmetric quantum mechanics (SUSYQM), initially proposed as a model for super- 

}y now developed into a field in its own 
8]. The essential ingredients of SUSYQM 
and charges Q- = (°[J), 



LtL o 
LLt 



symmetry (SUSY) breaking in field theory |7(, has 
right, with applications in diverse areas of physics 
are the supersymmetric Hamiltonian 7~t = y \ j 

Q + = fo^) which generate the supersymmetric algebra [H,Q±] = {Q±,Q±} = 0, 
{Q_,Q + } = Ti. The partner Hamiltonians Hi and H 2 satisfy an intertwining relation, 
LHi = H 2 L, where in one-dimension the transformation operator L = -jj + W(x) is a 
first-order Darboux transformation expressed in terms of a superpotential W(x). The inter- 
twining relation ensures that if is an eigenstate of H±, then also ^ 2 = L^>\ is an eigenstate 
of H 2 with the same energy, unless L^i vanishes or produces an unphysical state {e.g. non- 
normalizable) . Consequently, as shown in Fig. 1(a), the SUSY partner Hamiltonians Hi 
and H 2 are isospectral in the sense that their spectra consist of pair-wise degenerate levels 
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with a possible non-degenerate single state in one sector (when the supersymmetry is exact). 
The wave functions of the degenerate levels are simply related in terms of L and L'. Such 
characteristic features define a supersymmetric pattern. The intertwining relation ensures 
such properties for any pair of Hamiltonians not necessarily factorizable. We will continue 
to use the term "supersymmetric pattern" also in such circumstances. In what follows we 
focus the discussion on supersymmetric patterns obtained in selected Dirac Hamiltonians. 

The Dirac Hamiltonian, H, for a fermion of mass M moving in external scalar, Vs, and 
vector, Vj/, potentials is given by H = 6cp+j3(M +Vs) + Vv, where a, (3 are the usual Dirac 
matrices |9( and we have set K — c — 1. When the potentials are spherically symmetric: 
Vs = Vs(r), Vy = Vy(r), the operator K = —/3 (cr ■ £ + 1), (with er the Pauli matrices and 
I = —ir x V), commutes with H and its non-zero integer eigenvalues k = ±(j + 1/2) are 
used to label the Dirac wave functions ^ Kyin = ^ _1 (G' K [ Yex]m, iF K [Yf> x]m)- Here G K (r) 
and F K (r) are the radial wave functions of the upper and lower components respectively, 
and x are the spherical harmonic and spin function which are coupled to angular momentum 
j with projection m. The labels k = —(J + 1/2) < and £' = I + 1 hold for aligned spin 
j — £+1/2 (si/2,P3/2, etc.), while k = ( j + 1/2) > and £' = £ — 1 hold for unaligned spin 
j = £ — 1/2 (pi/2, g?3/2j etc.). Denoting the pair of radial wave functions by 

•■ = (£)■ (1) 

the radial Dirac equations can be cast in Hamiltonian form, H K $ K = E $ K , with 

( M + A -^- + - 

H K = dr r (2a) 

A(r) = Vs + V v , E(r) =V S - V v . (2b) 

We now look for Dirac Hamiltonians H Kl and H K2 which satisfy an intertwining relation of 
the form 

LH Kl = H K2 L . (3) 



Following 10( we consider a matricial Darboux transformation operator 



d 



L = A(r)— + B(r) 
dr 



(4) 



where A and B are 2x2 matrices with r-dependent entries Ay(r), B^{r). Relations (J3J) 
and (UJ) should be regarded as a system of equations for the unknown operator L and the 
so-far unspecified potentials in H K (j2J). The matrices A(r) and B(r) are found to be 



A 11 

E>22 

B21 



A 22 = a , A 12 
-b(M + A) - 



-A 
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1 „ 1 

— flW-K + 1) + -c 2 



b(M + S) - j-aw- (w+ - 1) + ^c 2 
aVV - ^buj + {uj_ + 1) + ^cx 



-aVV + — bu + (u- - r 
2r 



1 

2 Cl 



(5a) 
(5b) 
(5c) 
(5d) 
(5e) 



where uj + = (ki+Kz), ui- = («2 — Ki) and a, 6, ci, C2 are constants. In addition, the following 
relations have to be obeyed 



2a 



r J 



UJ-(UJ + + - 1) ^+ 







-41/ V (M + V S ) + 



w + (w_ + l)(w_ - 1) 



2b— [oj+(M + V s ) + V v 
r 



-ca 2ci (M + y 5 ) = 

r 



(6a) 



(6b) 



with denotes differentiation with respect to r. In the usual application of SUSYQM, one 
starts from a solvable Hamiltonian Hi and uses the intertwining relation to obtain a new 
solvable Hamiltonian H2. In the present case we employ a different strategy, namely, insist 
that both partner Hamiltonians H K1 and H K2 be of the form prescribed in Eq. (J2J with the 
same potentials, and look for solutions of Eq. (0) such that the potentials are independent of 
k. We find that there are six such solutions characterized by u;+,u;_ = 0, ±1. The solution 
with uj_ = is trivial (ki = k 2 ), L = —bH K + \c2l where I is the 2x2 unit matrix. The 
solutions with o>_ = ±1 lead to constant potentials Vg = Sq and Vy = Vq. The physically 



interesting solutions require u + = 0, 1, —1 and lead to the Coulomb, pseudospin and spin 
limits respectively. 

The Coulomb limit (ki + k 2 = 0) 
The solutions of Eq. (0) with uj + = K\ + n 2 = fix the potentials to be of Coulomb type 

Vs = y + So , V V = ^ + V , (7) 

with arbitrary strengths, as, ay. The constants So and Vq amount to constant shifts in the 
mass and Hamiltonian respectively, and henceforth will be omitted. In terms of the quanti- 



ties rji = {a s M + a v E) /\,r) 2 = {a s E + a v M) /A , A = VM 2 - E 2 , 7 = ^Jk 2 + a ; - n 



2 - o- 2 
s "y ' 



the quantization condition reads: 7 + ^1 = — n r (n r = 0, 1, 2, ... ), and leads to the eigen- 
values 9| Ef^] K /M = [—a s a v ± £y/ ^ 2 — of + a 2 / ]/(ay + ^ 2 ), where ^ = (n r + 7), and the 
K-dependence enters through the factor 7. The spectrum consists of two branches denoted 
by superscripts (+) and (— ). The corresponding eigenfunctions are 



where E = E^ K and F\ = F(—n r ,2'y + l,p), F 2 = F(—n r + 1, 27 + l,p) are confluent 
hypergeometric functions in the variable p = 2Ar. The states and energies in each branch 
are labeled by {n r ,K). It is also possible to express the results in terms of the principal 
quantum number N defined as N — n r + \k\, (N = 1, 2, . . . ). For n r > 1 the eigenvalues in 
each branch are two-fold degenerate with respect to the sign of k, i.e. E^ K = E^_ K and 
En r ] K = En r _ K . For n r = there is only one acceptable state for each k, which has k < for 
the (+) branch and k > for the (— ) branch. Equivalently, for a fixed principal quantum 
number N, the allowed values of Av £1X6 AC — ±1, ±2, . . . , ±(N - 1), -N for the (+) branch 
and k = ±1, ±2, . . . , ±(N — 1), +N for the (— ) branch of the spectrum. 

Focusing on the set of states with the levels are separated according 

to the value of \k\ = j + 1/2. For fixed k, Ef^] K is an increasing function of n r and, as 
shown in Fig. 1(b), for each value of j we have a characteristic supersymmetric pattern. 



There are two towers of energy levels, one for —\k\ (with n r = 0, 1, 2, ... ) and one for +\k\ 
(with n r = 1,2,.. .). The two towers are identical, except that the E^ =0 level at the 
bottom of the — \k\ tower is non-degenerate. Similar patterns of pair-wise degenerate levels 
with ±k appear also in the (— ) branch of the spectrum. However, since for fixed k, E^] k 
is a decreasing function of n r , the non-degenerate ^ ^ level is now at the top of the 
+ \k\ tower, resulting in an inverted supersymmetric pattern. From Eqs. ©-© we find the 
transformation operator to be 



_d i_ e+ _i_ Ma+ as_ _d , ay 

dr r fti Ki dr r 



\ 



as d ay_ _d £-_ _ Mat- 

ki dr r dr r Ki J 



(9) 



where e± = K\ + as<y±/ ' K\ and a± = (as ± ay)- The operator L connects degenerate states 
with (n r > 1, ±ac), and annihilates the non-degenerate states with n r = 

L$ nr;Kl = C$ nr)K2 (ki = -k 2 ) ■ (10) 



Here C = ^ynAj ~ Vi) an d $n r ,K are given in Eq. (JHJ). Constructing supersymmetric 
charges Q± and Hamiltonian 7i from L and -ff Kl , -ff K2 in the manner described at the 
beginning of the letter, ensures that [T~t,Q±] = {Q±,Q±} = 0, but now {Q-,Q+} oc 
(H — E^ =0f .)(H — E^ =0k ), with k = Ki = —k 2 . These relations represent a quadratic 



3, 



which arises because both the 



deformation of the conventional supersymmetric algebra 
Dirac Hamiltonian and the transformation operator L are of first order. 

The explicit solvability and observed degeneracies of the relativistic Coulomb problem 
are related to the existence of an additional conserved Hermitian operator 

B = -iK l5 (H - /3M) + ^—^ ( a v M + a s H ) , (11) 

which commutes with the full Dirac scalar and vector Coulomb Hamiltonian, H, but an- 
ticommutes with K. This operator is a generalization of the Johnson-Lippmann operator 
applicable for as = [ill ]. 
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The pseudospin limit (k% + k 2 = 1) 
The solutions of Eq. © with uj + = K\ + k 2 = 1 require that the sum of scalar and vector 
potentials is a constant 



A(r) =V s (r)+V v (r) = A . 



(12) 



Under such condition the Dirac Hamiltonian is invariant under an SU(2) algebra, whose 



generators are 
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13] 



S, 



Sp 



s, 



(13) 



/ 



Here Sp = cr fJi /2 are the usual spin operators, Sp = U p SpU p and U p = — . This relativistic 



symmetry has been used |3( to explain the occurrence of pseudospin doublets in nuclei 
The latter refer to the empirical observation of quasi-degenerate pairs of certain normal- 
parity shell-model orbitals with non-relativistic single-nucleon radial, orbital, and total an- 
gular momentum quantum numbers: (n, £,j = i + 1/2) and (n — 1,£ + 2, j — £ + 3/2). The 
doublet structure is expressed in terms of a "pseudo" orbital angular momentum, I — £+ 1, 
and "pseudo" spin, s = 1/2, which are coupled to j = £±s. For example, (nsi/2, (n — 1) dx/o) 
will have I — 1, etc. Such doublets play a central role in explaining features of nuclei |l5|. 
including superdeformation and identical bands. In a relativistic description of nuclei |l|, 
these non-relativistic wave functions are identified with the upper components of Dirac wave 
functions, ^ Kl <o,m an d ^K 2 >o,m with K\ + k 2 = 1, which are eigenstates of a Dirac Hamilto- 
nian with scalar and vector mean field potentials, approximately satisfying condition (|12|). 



The corresponding lower components have n nodes [16| and orbital angular momentum 
equal to I j3j. In the pseudospin limit these two Dirac states form a degenerate doublet, 



and their radial components satisfy F K1 = F K2 , and —p 1 - 



dG^2 
dr 



+ ^G K2 . These 



relations have been tested in numerous realistic mean field calculations in a variety of nu- 
clei, and were found to be obeyed to a good approximation, especially for doublets near the 



Fermi surface 
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Q. 



For potentials with asymptotic behavior as encountered in nuclei, 



the Dirac eigenstates for which both the upper (G K ) and lower (F K ) components have no 
nodes, can occur only for k < 0, and hence do not have a degenerate partner eigenstate 



(with k > 0) [16[. These nodeless Dirac states correspond to the shell-model states with 
(n = 0,£,j = £ + 1/2). For heavy nuclei such states with large j, i.e., 0g 9 / 2 , 0h u / 2 , 0zi 3 / 2 , 
are the "intruder" abnormal-parity states which, indeed, empirically are found not to be 



part of a doublet |15(. Altogether, as shown in Fig. 1(c), the ensemble of Dirac states with 
«4 + k>2 — 1 exhibits a supersymmetric pattern of twin towers with pair-wise degenerate 
pseudospin doublets sharing a common £, and an additional non-degenerate nodeless state 
at the bottom of the n\ < tower. An exception to this rule is the tower with k 2 = 1 (pi/2 
states with £ = 0), which is on its own, because states with K\ = do not exist. From 
Eq. ©-(0) we find the transformation operator to be 

( n d Ko \ 

(14) 







_d K2. 

dr r 



_d_ 
\ dr 



J 



f (2M + S + A ) 

L connects the two doublet states as in Eq. (fTUj) but with ki + k 2 = 1 and C = b(M+A — E). 



In this case, {Q-, Q + } = b 2 [H — (M + A )][7Y — (Af + A )] is proportional to a polynomial 
of 7i, again indicating a deformation of the conventional SUSY algebra. In real nuclei, the 
relativistic pseudospin symmetry is slightly broken, implying A(r) ^ Ao in Eq. (fT2"|) . Taking 
H K as in Eq. (0) and L as in Eq. (fT4^) but with A — > A(r), we find that LH K1 — H K2 L = 
ibA'<j2- Furthermore, {Q-, Q+} has the same formal form as before, but the appearance of 
A(r) instead of A implies that the anticommutator is no longer just a polynomial of 7i. 

The spin limit (ki + k 2 = — 1) 
The solutions of Eq. @ with uj + = n\ + k 2 = — 1 require that the difference of the scalar 
and vector potentials is a constant 



Ur) = Vs(r) - Vv(r) = E n . 



(15) 



Under such condition the Dirac Hamiltonian is invariant under another SU(2) algebra, whose 



S 



generators are obtained from Eq. (jlHj) by interchanging s M and s M 



s M 



(16) 



This relativistic symmetry gives rise to degenerate doublets of Dirac states ^ Kl< o, m and 
^«;2>o,m with Ki + k 2 = — 1, whose upper components have quantum numbers (n,£,j = 
£+1/2) and (n, £, j = £—1/2). Such spin doublets were argued to be relevant for degeneracies 
observed in heavy-light quark mesons |4] and possibly for the anti-nucleon spectrum in a 
nucleus 1^. In the spin limit, the corresponding radial components satisfy G Kl = G K2 and 
— ^F Kl = -4^ — ^F K2 . As shown in Fig. 1(d), the spectrum consists of towers of states 
with Ki + k 2 - 



-1, forming pair-wise degenerate spin doublets. In this case, none of the 
towers have a single non-degenerate state and hence, the spectrum corresponds to that of 
a broken SUSY -8 1 ]. The tower with K\ = —1 (si/2 states) is on its own, since states with 
k 2 = do not exist. The transformation operator, found from Eqs. 



/ 



V 



(2M + S + A) 



d I Ki ^ 

dr r 



_a I K2 

dr r 







(17) 



connects the two doublet states as in Eq. (fTU|) . but with K\ + k 2 = —1 and C = —b(E + 
M + S ). The nilpotent charges, Q±, commute with the supersymmetric Hamiltonian, 7i, 
and exhibit a deformation of the conventional SUSY algebra, {Q_,Q + } = b 2 [H + (M + 
E )][W + (M + E )]. When S(r) + S in Eq. (USD, we have LH K1 — H K2 L = -ibE'a 2 , where 
H K is given in Eq. (J2J) and L, as well as {Q-, Q+}, have the same form as before but with 
S - S(r). 

In summary, a common intertwining relation was shown to provide the basis for a unified 
treatment of three separate limits at which a Dirac Hamiltonian, with scalar and vector 
potentials, exhibits supersymmetric patterns. In the Coulomb limit the potentials are 1/r 
but their strengths are otherwise arbitrary. In the pseudospin or spin limits there are no 
restrictions on the r-dependence of the potentials but there is a constraint on their sum 
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131 ] . It is gratifying to note that 



or difference. The characteristic degeneracies reflect the presence of additional conserved 
operators, the generalized Johnson-Lippmann operator given in Eq. fill}) , and the previously 
introduced relativistic pseudospin and spin generators 
the indicated supersymmetric patterns are manifested empirically, to a good approximation, 
in physical dynamical systems. While previous studies have focused on individual doublets 
in nuclei and hadrons, it is the grouping of several doublets (and intruder levels in nuclei), as 
suggested in the present work, which highlights the fingerprints of supersymmetry present 
in these dynamical systems. This work was supported by the Israel Science Foundation. 



10 



[1] B.D. Serot and J.D. Walecka, Adv. Nucl. Phys. 16, 1 (1986). 

[2] W. Lucha, F.F. SchSberl and D. Gromes, Phys. Rep. 200, 127 (1991). 

[3] J.N. Ginocchio, Phys. Rev. Lett. 78, 436 (1997). 

[4] P.R. Page, T. Goldman and J.N. Ginocchio, Phys. Rev. Lett. 86, 204 (2001). 
[5] B. Thaller, The Dirac Equation, (Springer- Verlag, 1992). 
[6] C. V. Sukumar, J. Phys. A 18, L697 (1985). 
[7] E. Witten, Nucl. Phys. B 188, 513 (1981). 

[8] G. Junker, Supersymmetric Methods in Quantum and Statistical Physics, (Springer- Verlag, 
Heidelberg 1996). 

[9] W. Greiner, B. Miiller and J. Rafelski, Quantum Electrodynamics of Strong Fields (Springer- 
Verlag, NY, 1985). 

[10] N. Debergh et al, J. Phys. A 35, 3279 (2002); L.M. Nieto et al, Ann. Phys. 305, 151 (2003). 

[11] M. H. Johnson and B. A. Lippmann, Phys. Rev. 78, 329 (1950). 

[12] J.S. Bell and H. Ruegg, Nucl. Phys. B 98, 151 (1975). 

[13] J.N. Ginocchio and A. Leviatan, Phys. Lett. B 425, 1 (1998). 

[14] K.T. Hecht and A. Adler, Nucl. Phys. A 137, 129 (1969); A. Arima, M. Harvey and K. 

Shimizu, Phys. Lett. B 30, 517 (1969). 
[15] A. Bohr, I. Hamamoto and B.R. Mottelson, Phys. Scr. 26, 267 (1982). 
[16] A. Leviatan and J.N. Ginocchio, Phys. Lett. B 518, 214 (2001). 

[17] J.N. Ginocchio and D.G. Madland, Phys. Rev. C 57, 1167 (1998); G.A. Lalazissis et al, ibid. 

58, R45 (1998); J. Meng et al, ibid. 59, 154 (1999). 
[18] J.N. Ginocchio and A. Leviatan, Phys. Rev. Lett. 87, 072502 (2001); J.N. Ginocchio, Phys. 

Rev. C 66, 064312 (2002). 
[19] S.G. Zhou, J. Meng and P. Ring, Phys. Rev. Lett. 91, 262501 (2003). 

11 



FIG. 1: Schematic qualitative supersymmetric patterns in (a) SUSYQM and in the (b) Coulomb, 
(c) pseudospin, (d) spin, limits of the Dirac Hamiltonian. In (a) H\ and H2 have identical spectra 
with an additional level for H\ when SUSY is exact. Spectroscopic notation nij in (b)-(d) refers 
to quantum numbers of the upper component, and re, N, i are defined in the text. In (b) the radial 
nodes n are related to n r by n r = n (n r = n + 1) for k < (re > 0), and only the E^] K branch is 
shown. 
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